Crystal resonant frequency sensor

ABSTRACT

A method and apparatus for measuring the frequency of a desired resonant mode of a crystal arrangement, or other two-port device, during an automated operation. The crystal arrangement, or other two-port device, is placed into a test circuit and subjected to a sinusoidal test signal of known frequency. Based upon the output response of the crystal arrangement to the test signal, the frequency of the test signal is changed such that the test signal rapidly converges on a desired mode of operation of the crystal arrangement. This is accomplished by first noting a desired increase in amplitude of the output response of the crystal arrangement, followed by measuring an error signal related to the desired crystal arrangement mode of operation. When the error equals a predetermined value the frequency of the sinusoidal test signal is the frequency of the desired mode.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to frequency measurement and, more particularly, to frequency measurements of various quartz crystal resonant modes.

2. Description of the Prior Art

It is often desirable to determine the resonant modes of a quartz crystal arrangement where the frequency information may be used in scientific study, for example. Frequency characteristics of interest may include resonant frequencies, inharmonic overtones, and harmonic modes of the crystal arrangement.

Historically, methods of measuring resonant modes of a crystal include swept frequency sinusoidal analysis and crystal oscillator testing. The swept frequency sinusoidal analysis involves applying a sinusoidal signal to the crystal arrangement and sweeping the frequency of the sinusoidal input signal through a frequency range of interest. Based on the response of the crystal arrangement to the sinusoidal input, certain frequency characteristics can be obtained. The swept frequency sinusoidal analysis, however, is often tedious and overly time consuming, and often requires repetitive cycles of: providing an input signal to the crystal arrangement, changing the frequency of the input signal, measurement of the response of the crystal arrangement, comparison of the gathered information, and determining the input signal to be applied to the crystal arrangement for a subsequent measurement based on the current cycle measurement. Typically, for a high accuracy measurement very precise and expensive measurement equipment must be used. However, even when measurements are taken with the very precise and expensive equipment, errors may still arise due to human intervention during the analysis.

The crystal oscillator testing method for measuring a resonant mode of a crystal involves placing the crystal in an electronic oscillator circuit and measuring the oscillator signal frequency. Due to the resonant properties of a crystal, the electronic oscillator circuit will converge to the crystal's resonant frequency. While crystal oscillator testing is a reasonable method to determine the principal resonant frequency of a crystal, it is an undesirable method where additional resonant modes are to be measured.

OBJECTS OF THE INVENTION

It is one object of this invention to provide a method and apparatus to rapidly and accurately determine the frequency of a desired resonant mode characteristic of a crystal arrangement, or other two port device, under test without the use of complicated and expensive test equipment.

It is another object of this invention to provide a method and apparatus to rapidly determine the frequency of a desired resonant mode of a crystal arrangement, or other two port device, under test as part of a fully automated operation.

It is another object of this invention to provide a method and apparatus to rapidly determine the frequency of a desired resonant mode of a crystal arrangement under test while discriminating between several different modes of oscillation of the crystal.

It is another object of this invention to provide a method and apparatus for measuring multiple, closely spaced resonant frequencies of a crystal arrangement, or other two port device, under test.

It is another object of this invention to provide a precision, linear crystal-controlled oscillator utilizing an inexpensive crystal and a voltage controlled oscillator.

SUMMARY OF THE INVENTION

The foregoing and other problems are overcome and the objects of the invention are realized by methods and apparatus in accordance with embodiments of this invention. In accordance with one embodiment of the invention an apparatus and method are provided for measuring desired resonant modes of a crystal arrangement, or other two-port device. First, based upon calculated operating parameters related to a desired resonant mode of the crystal arrangement, the apparatus and method measure the passive response of the crystal arrangement to a sinusoidal input of a known frequency. The passive response of the crystal arrangement is measured by sweeping the frequency until the amplitude of the response is above a calculated threshold value. Second, the apparatus and method measure an error signal. The error signal represents the difference in frequency between the sinusoidal input and the desired resonant mode, generated, in part, as a result of the desired resonant mode of the crystal arrangement. Additionally, the apparatus and method fine tune, or closed-loop feedback adjust the frequency of the sinusoidal input such that the error signal rapidly converges to a predetermined value. When the error signal reaches the predetermined value, the frequency of the sinusoidal input is the frequency of the desired resonant mode of the crystal arrangement.

In accordance with the present invention, a crystal resonant frequency sensor includes a controller whose output is responsive to values measured at a minimum of two inputs, a voltage-controlled oscillator adapted to provide a signal whose frequency is responsive to a command signal, a phase shifter, a first and second multiplier and a first and second filter.

BRIEF DESCRIPTION OF THE DRAWINGS

The above set forth and other features of the invention are made more apparent in the ensuing Detailed Description of the Invention when read in conjunction with the attached Drawings, wherein:

FIG. 1 is a schematic diagram of an embodiment of a crystal resonant frequency sensor in accordance with this invention.

FIG. 2 is a graph that illustrates two signals of a crystal resonant frequency sensor, plotted with respect to normalized frequency.

FIG. 3 is a flow chart depicting a first part of a method in accordance with this invention.

FIG. 4 is a flow chart depicting a second part of a method in accordance with this invention.

FIG. 5 is a flow chart depicting a third part of a method in accordance with this invention.

FIG. 6 is a graph that illustrates a first signal, the amplitude squared of oscillation, of a crystal resonant frequency sensor, plotted with respect to time.

FIG. 7 is a graph that illustrates a second signal, the error or detuning in frequency, of a crystal resonant frequency sensor, plotted with respect to time.

FIG. 8 is a graph that illustrates a normalized signal, h of a crystal circuit, plotted with respect to frequency.

FIG. 9 is a schematic diagram of a two-port mechanical structure whose resonant characteristics can be measured by the present invention.

FIG. 10 is a graph that illustrates a normalized signal, h of a two-port mechanical system, plotted with respect to frequency around the first resonant peak.

FIG. 11 is a graph that illustrates a normalized signal, g of a two-port mechanical system, plotted with respect to frequency around the second resonant peak.

FIG. 12 is a graph that illustrates a normalized signal, g of a two-port mechanical system, plotted with respect to frequency around the second resonant peak.

DETAILED DESCRIPTION OF THE INVENTION

Referring to the accompanying drawings, FIG. 1 depicts an embodiment of a crystal resonant frequency sensor. Specifically, FIG. 1 depicts a crystal resonant frequency sensor 10 comprising a controller 14, a voltage-controlled oscillator (VCO) 16, a phase shifter 18, a first multiplier 20, a second multiplier 22, a first filter 24, and a second filter 26. A crystal arrangement 12 to be characterized is shown connected between the output 12 a of the phase shifter 18 and a ground 12 b. The voltage e₂ is measured with respect to ground at the output 12 a of the phase shifter 18. The crystal arrangement 12 has an input 12 a and an output 12 b. The crystal arrangement 12 may comprise a crystal element which inherently comprises stray capacitance (C_(d)), series capacitance (C), a value Q representing a relationship between energy stored and energy dissipated per cycle, and at least one resonant mode. The crystal arrangement 12 may be modeled as a one mode equivalent circuit, as shown, comprising stray capacitance (C_(d)) coupled in parallel to a series combination of inductance (L), parallel capacitance (C) and resistance (r), or the crystal arrangement 12 may generally comprise components forming a circuit which has at least one resonant mode.

The controller 14 further comprises a first input 14 a, a second input 14 b, and an output 14 c. Optionally, the controller 14 may further include an input 14 d for user inputted parameters. The first multiplier 20 and the second multiplier 22 both further comprise a first input 20 a, 22 a, and a second input 20 b, 22 b, and an output 20 c, 22 c, respectively. The output 14 c of the controller 14 electrically communicates with a control signal input of the VCO 16. The output e₁ of the VCO 16 is electrically coupled with an input to the phase shifter 18, and the first input 22 a of the second multiplier 22. The output 12 a of phase shifter 18 is electrically coupled to the second input 22 b of the second multiplier 22, the first input 20 a of the first multiplier 20, and the second input 20 b of the first multiplier 20. The crystal arrangement 12 under test is electrically connected in shunt with the phase shifter 18, between the phase shifter 18 and the multipliers 20, 22 as shown. Thus, the signal e₂ is generated, in part, by the passive response of crystal arrangement 12 to the output signal e₁ of the VCO 16.

The first filter 24 and the second filter 26 both comprise an input and an output. The output 20 c of the first multiplier 20 is electrically coupled to the input of the first filter 24. The output 22 c of the second multiplier 22 is electrically coupled to the input of the second filter 26. Additionally, the output of the first filter 24 is electrically coupled to the input 14 a of the controller 14. The output of the second filter 26 is electrically coupled to the input 14 b of the controller 14. The input 14 a of the controller 14, denoted h in FIG. 1, represents a signal which is proportional to the amplitude or modulus of the signal e₂. The input 14 b of the controller 14, denoted g in FIG. 1, represents an error directly related to the phase relationship of the signal e₁ with respect to the signal e₂.

In operation, the sensor 10 operates in either an open-loop mode or a closed-loop mode. The sensor 10, after an initialization sequence, operates in the open-loop mode. As a precursor to initiating the open-loop mode of operation a user inputs one or more parameters related to a desired resonant mode of the crystal arrangement 12 under test. The inputted parameters may include values of Q, a ratio of capacitance representing a relationship between stray capacitance (C_(d)) across the crystal and crystal series capacitance (C), and the approximate value of the resonant frequency related to the desired resonant mode. The user may input the parameters using any suitable means, such as keypad 14 d that forms a part of controller 14, or a computer (not shown) in serial communication with a data port of the controller 14. Based on the user input, the controller 14 calculates operating parameters which may comprise a gain (K) which the controller 14 applies to the input value of h and g, a threshold value (h_(th)) which is to compared to the input value of h, and sweep parameters, including initial and peak oscillator voltage values corresponding to the desired frequency range of operation along with a starting point frequency and a frequency sweep rate. It has been discovered through testing that only approximate values of the above parameters are needed in order to acquire the frequency of the desired mode of a crystal arrangement 12 under test. However, the ability of the crystal resonant frequency sensor to rapidly converge on the frequency of a desired mode is directly related to the calculated parameters, as will become apparent.

The controller 14, the VCO 16, the phase shifter 18, the first multiplier 20, the first filter 24, and the crystal arrangement 12 comprise a control loop during the open-loop mode of operation. In the open-loop mode the input 14 a (h) of the controller 14 is enabled and input 14 b (g) of the controller 14 is ignored. Generally, in the open-loop mode a sinusoidal signal e₁ with a frequency within the calculated frequency range is provided, by way of the phase shifter 18, to the crystal arrangement 12. In response to this sinusoidal input, the crystal arrangement 12 passively provides an output signal e₂. The first multiplier 20 and the first filter 24 generate a signal h which is proportional to crystal arrangement 12 output signal e₂. The frequency of the sinusoidal signal e₁ is then adjusted until the magnitude of the value of h exceeds the calculated threshold value h_(th). An increasing value of h is an indication that the sinusoidal frequency of signal eis approaching the desired resonant frequency of the crystal arrangement 12. Thus, when the value of h exceeds the predetermined threshold value h_(th) the corresponding frequency of signal e₁ is assumed to be approaching the desired resonant frequency of the crystal arrangement 12, and the controller 14 initiates the closed-loop mode of operation.

More specifically, the controller 14 initiates the open-loop mode by providing a signal 14 c to the VCO 16 commanding the VCO 16 to output a sinusoidal signal e₁ at a frequency which is within the calculated frequency range such that the value of h initially does not exceed the calculated threshold value h_(th). Typically, the initial frequency of sinusoidal signal e₁ is at a limit of the calculated frequency range, the lower frequency limit of the range, for example. While the phase shifter 18 is physically in the path of the open-loop mode of operation, it is of little concern during the open-loop mode of operation since the amplitude of the output signal e₂ of the phase shifter 18 is unchanged by the phase shifter 18 itself. However, in response to the signal e₁, the crystal arrangement 12 passively provides the signal e₂. Thus, signal e₂ is related to the frequency characteristics of the crystal arrangement 12.

The first multiplier 20 and the first filter 24 work together to act as an amplitude demodulator for the signal e₂. One skilled in the art of electronics will recognize that the multiplier produces signals near zero frequency and at twice the input frequency that are proportional to the amplitude of the input signal. The low-pass filter allows a direct current (DC) representation of the signal to be available for the further operations. In a preferred embodiment, the first filter 14 is a low pass filter. The low pass filter has a cutoff frequency for h set at a fraction of the resonant frequency (ω_(o)). Note, the specific design of the sensor 10 determines the cutoff frequency of the low pass filter, which must exclude ω_(o) and 2ω_(o), and must be large enough to pass the frequency sweep rate transients. The controller 14 then measures the signal h and compares the measured value of h with the predetermined threshold value h_(th). If the measured value of h is less than the threshold value h_(th), the controller 14 adjusts the output 14 c to command the VCO 16 to adjust the frequency of the signal e₁ such that the value of h approaches the threshold value h_(th).

Referring momentarily to FIG. 2, a typical normalized steady-state response of h versus frequency is graphically shown. For clarity the curve representing h is shown as a dashed line. As described above, the initial frequency selection for signal e₁ is calculated to produce an initial value of h which is below the calculated threshold value h_(th), point A on the graph, for example. As the controller 14 sweeps the frequency of signal e₁ towards the desired resonant frequency of the crystal arrangement 12 under test, the value of h will approach the value of h_(th), depicted by point B, in the direction of arrow h₁. When the value of h surpasses the value of h_(th), depicted by point C, the crystal resonant frequency sensor enters the closed-loop mode of operation. While this example depicts the response of the present invention to an initial value of h corresponding to a frequency which is lower than the frequency of the desired mode, the outcome would be similar for an initial value of h whose corresponding frequency is greater than the frequency of the desired mode. As discussed above, the controller 14 commands the VCO 16 to adjust the frequency of the signal e₁ such that the value of h approaches the threshold value h_(th). If the initial value of h corresponded to a frequency which was greater than the frequency of the desired mode of crystal arrangement 12, shown as point A1, the controller 14 would then command the VCO 16 to decrease the frequency of signal e₁. The decrease in frequency of signal e₁ causes the value of h to move towards the value of h_(th) along the path indicated by arrow h₂. As discussed above, when the value of h surpasses the value of h_(th) the crystal resonant frequency sensor enters the closed-loop mode of operation.

It is important to note that the controlled operation of the VCO 16 can be utilized in preliminarily determining parameters of different crystal arrangements. For example, for particular crystal arrangements where the crystals are in contact with a liquid or suspended in a liquid, the oscillator can be commanded to sweep through a frequency range and the corresponding value of h can be measured. This information can be used during a subsequent measurement of a resonant mode frequency by the crystal resonant frequency sensor. The negative feedback nature of the crystal resonant frequency sensor is an advantage since it is less sensitive to varying parameters of crystal arrangements, such as the one described above, than a positive feedback system, such as an oscillator.

Turning back to FIG. 1, the controller 14, the VCO 16, the phase shifter 18, the second multiplier 22, the second filter 26, and the crystal arrangement 12 comprise the control loop during the closed-loop mode of operation. In the closed-loop mode the controller 14 adjusts the frequency of the VCO 16 in order to drive g towards zero. As the measured value of g approaches the predetermined value, the corresponding frequency of the VCO 16 approaches the resonant frequency of the crystal arrangement 12.

More specifically, as stated above, once the value of h as measured by the controller 14 is greater than the predetermined value h_(th), the closed-loop mode of operation is initiated. During the closed-loop mode of operation the controller also enables input 14 b (g). This signal (g) is the primary control signal, but 14 a (h) is used during the closed-loop mode to test for convergence. When the frequency of signal e₁ approaches the resonant frequency of the crystal arrangement 12, the phase of signal e₂ with respect to signal e₁ approaches 90° and the resultant output 20 c (g) of the second multiplier approaches a zero value. The output 22 c of the second multiplier 22 is then provided to the input of the second filter 26 where unwanted frequency components of the signal are filtered. In a preferred embodiment the second filter 26 is a low pass filter whose cutoff frequency for g must pass the same requirements as h but additionally be significantly larger than the bandwidth of the servo loop. The output of the second filter 26 is provided as input 14 b (g) to the controller 14. When the absolute value of the input signal 14 b (g) of controller 14 is measured to be less than some convergence value, the corresponding frequency of signal e₁ is the desired resonant frequency of the crystal arrangement 12 under test. The resonant frequency can be determined, for example, from the voltage used for controlling the VCO, or from measurements of the output of the VCO using a counter. For low accuracy applications, the resonant frequency can be determined from the VCO control voltage, after calibration of the VCO. For high accuracy applications, one of two methods for determining resonant frequency can be utilized. First, an expensive VCO is used where the VCO voltage is calibrated and the resonant frequency measurement is derived from the VCO voltage. Second, an ordinary VCO is used with a precise counter which can determine the resonant frequency.

Turning momentarily back to FIG. 2, the transitional phase between the open-loop mode and the closed-loop mode of operation can be more easily understood. When the value of h at the input of 14 a of controller 14 exceeds the value h_(th), denoted as point C of the graph, the controller 14 enters the closed-loop mode and ignores input 14 a (h) while enabling input 14 b (g). It should be noted that while the point C is shown clearly above the value of h_(th), this is done for purposes of clarity. The actual location of point C is dependent on the resolution of the acquisition of h by the controller 14, defined in part by the gain (K) of the system. The controller 14 measures the initial value of input 14 b (g) at this time, denoted g_(int) at point D of the graph. As shown, the signal g is an approximately linear function over the crystal resonant bandwidth. The controller 14 then finely adjusts the frequency of signal e₁ such that the measured value of g at input 14 b of controller 14 rapidly approaches a zero value, denoted point E of the graph. As discussed above, when the value of g is measured to be approximately zero, the frequency of the signal e₁ represents the desired resonant frequency of the crystal arrangement 12 under test.

If, however, the initial value of h corresponds to a signal e₁ frequency greater than the desired resonant frequency, the initial value of g will be positive. The controller 14, attempting to minimize g, decreases the frequency of signal e₁.

At this point it is important to note that the phase shifter 18 may be replaced by a generalized network in order to determine the desired resonant modes of other two-port resonant devices. Like the crystal arrangement discussed above, such two-port devices would comprise an input which receives the e₁ signal from the generalized network, and an output which would be measured to determine the signal e₂. The loop dynamics would be consistent with that stated above wherein the signal g represents an error function of the frequency error between the resonant frequency and the oscillator frequency with the value equaling a predetermined value, such as a zero value, at zero error. The generalized network would be designed such that the error function converges to the desired predetermined value.

Additionally, it would be apparent to one of ordinary skill in the art to combine all or a part of the elements described above into a single element. For example, microcontrollers or other similar microprocessor based devices have become more readily available in the past decade. Such devices, for example a controller chip from the Motorola M68HCXX family of products, typically have standard features that provide I/O for receiving digital and analog information and memory for user programs and standard functions. The selection of a particular controller chip, for example an eight bit chip (M68HC11) versus a sixteen bit chip (M68HC16), depends on the accuracy required for the frequency command to the VCO 16. Once selected, the controller chip could be programmed to provide a command to the VCO 16 to excite the crystal with a sinusoidal signal of a given frequency. As such, a suitable digital signal processor (DSP) device can be used to implement a number of the discrete functions shown in FIG. 1.

Assuming a purely capacitive phase shifter with capacitance, C₀, the crystal arrangement 12, which as shown is an equivalent representation of an actual crystal device, may be modeled by the following normalized first order equations which describe the four element resonant circuit of the crystal arrangement 12: $\begin{matrix} {\frac{v}{t} = {{{- \omega_{o}^{2}}q} - {\frac{\omega_{o}}{Q}v} + {\mu \quad e_{1}}}} & (1) \\ {\frac{q}{t} = v} & (2) \\ {e_{2} = \frac{\left( {{C_{O}e_{1}} - q} \right)}{\left( {C_{O} + C_{d}} \right)}} & (3) \end{matrix}$

Where:

C_(O)=the capacitance across phase shifter 18; $\begin{matrix} {{\omega_{O}^{2} = \frac{\left( {C_{O} + C_{d} + C} \right)}{{LC}\left( {C_{O} + C_{d}} \right)}};} \\ {{\mu = \frac{C_{O}}{L\left( {C_{d} + C_{O}} \right)}};{and}} \\ {Q = {\frac{\left( {L/C} \right)^{1/2}}{r}\left\lbrack \frac{C_{O} + C_{d} + C}{C_{O} + C_{d}} \right\rbrack}^{1/2}} \end{matrix}$

The variables q, v, e₁, and ω_(o) may be further defined as follows: $\begin{matrix} \begin{matrix} {q = \quad {C\left\lbrack {{{\xi cos\omega}\quad t} + {\eta \quad \sin \quad \omega \quad t}} \right\rbrack}} \\ {v = \quad {C\quad {\omega \left\lbrack {{{- \xi}\quad \sin \quad \omega \quad t} + {\eta \quad \cos \quad \omega \quad t}} \right\rbrack}}} \\ {e_{1} = \quad {E_{1}\cos \quad \left( {{\omega \quad t} + \phi} \right)}} \\ {\omega_{O} = \quad {\omega + {\Delta \quad \omega_{O}}}} \end{matrix} & (4) \end{matrix}$

where the variables ξ and η describe a response envelope of the crystal arrangement 12. Substituting the defining equations (4) for q and v into the differential equations of (1) and (2) above results in a pair of differential equations in terms of the envelope parameters ξ and η. The envelope variables ξ and η are then assumed to vary slowly compared to a period of oscillation and the following differential equations for the envelope are obtained: $\begin{matrix} {\begin{matrix} {\frac{\xi}{t} = {{\Delta \quad \omega_{O}\eta} - {\frac{\omega_{O}}{2Q}\xi} + {\frac{\rho \quad E_{1}}{2}\sin \quad \phi}}} \\ {\frac{\eta}{t} = {{{- \Delta}\quad \omega_{O}\xi} - {\frac{\omega_{O}}{2Q}\eta} + {\frac{\rho \quad E_{1}}{2}\cos \quad \phi}}} \end{matrix}{{{Where}\text{:}\quad \rho} = \frac{C_{O}}{\left( {C_{O} + C_{d} + C} \right)}}} & (5) \end{matrix}$

The error signal g and the signal h, can be characterized in terms of the response of the crystal arrangement 12 to the voltage-controlled oscillator 16 excitation by the following: $\begin{matrix} \begin{matrix} {g = {\frac{E_{1}^{2}}{2} - {\left\lbrack \frac{C}{2C_{O}} \right\rbrack {E_{1}\left( {{\xi \quad \cos \quad \phi} - {\eta \quad \sin \quad \phi}} \right)}}}} \\ {h = {\frac{E_{1}^{2}}{2} + {{\frac{1}{2}\left\lbrack \frac{C}{C_{O}} \right\rbrack}^{2}\left( {\xi^{2} + \eta^{2}} \right)} - {\left\lbrack \frac{C}{C_{O}} \right\rbrack {E_{1}\left( {{\xi \quad \cos \quad \phi} - {\eta \quad \sin \quad \phi}} \right)}}}} \end{matrix} & \text{(6)} \end{matrix}$

The control dynamics of the voltage-controlled oscillator 16 can be modeled as a sinusoidal signal with a given phase, as follows: $\begin{matrix} {e_{1} = {E_{1}\cos \quad \left( {{\omega \quad t} + \phi} \right)}} & (7) \end{matrix}$

Where:

ω is a constant “nominal frequency” parameter; and

e₁ is the voltage-controlled oscillator output voltage.

Further, the change of phase, since the phase is a slowly varying function of time as compared with the period of oscillation, can be represented by: $\begin{matrix} {\frac{\phi}{t} = {\Omega + \Omega_{n}}} & (8) \end{matrix}$

Where:

Ω is the command to the voltage-controlled oscillator in radian frequency; and

Ω_(n) the frequency noise referred to the input of the voltage-controlled oscillator 16.

The feedback control law corresponding to the closed-loop mode of operation around the error signal g is defined as: $\begin{matrix} {\frac{\Omega}{t} = {- {Kg}}} & (9) \end{matrix}$

Where K is the gain of the feedback loop.

Equation (9) provides an integration so that a constant error in frequency g is integrated with respect to time and thus increases the command Ω to the voltage-controlled oscillator 16 until feedback decreases the error g and the command Ω settles to a constant value. The command Ω represents a proportional frequency change in the output signal e₁ of the voltage-controlled oscillator 16.

In order to determine how close the voltage-controlled oscillator 16 output frequency is to the resonant frequency of the crystal arrangement 12 the steady state response of the crystal resonant frequency sensor 10 is calculated from equations (5) through (9). In steady-state the phase Ω advances with constant frequency Ω_(s) and the value of the error signal g is equal to zero. Using classical sinusoidal methods, equations (5) through (9) above are solved to obtain the steady state condition for the error signal g: $\begin{matrix} {g = {{\frac{E_{1}^{2}}{2} - {\left( \frac{C}{C_{O}} \right){\left( \frac{\rho \quad \omega_{O}E_{1}^{2}}{2} \right)\left\lbrack \frac{\Omega_{s}}{\Omega_{s}^{2} + \left( {{\omega_{O}/2}Q} \right)^{2}} \right\rbrack}}} = 0}} & (10) \end{matrix}$

Equation (10) then must be solved for Ω_(s) in order to determine the corresponding resonant frequency of the crystal arrangement 12.

Referring now to FIG. 3, FIG. 4, and FIG. 5, a flow chart 30 is shown depicting a method in accordance with this invention. Specifically, FIG. 3 shows an initialization part of the method depicted in the flow chart 30 initiated with a first start step 32. Step 32 is followed by a step 34 where the user enters one or more parameters, as described above, into the controller 14 via an input device such as the keypad 14 d. Upon receiving the user input the controller 14 calculates the operating parameters in a step 36. Once the operating parameters are calculated the open-loop mode of operation is initiated in a step 38.

Following the circled-1 to FIG. 4, in a step 40 the controller 14 commands the VCO 16 to provide the signal e₁ at a given frequency which was determined in step 36. The value of the amplitude signal h is then measured by the controller 14 in a step 42. In a step 44 the controller 14 compares the value of h with a predetermined value h_(th) calculated in the step 36. If the value of h is less than the value of the signal h_(th) the controller 14 commands the VCO 16 to adjust the frequency of signal e₁ in a step 46. Steps 42, 44, and 46 are repeated until the value of the signal h is greater than the value of the signal h_(th). If, during the step 44, the value of signal h is determined to be greater than the value of the signal h_(th) by the controller 14, the controller 14 initiates the closed-loop mode of operation in a step 48.

Following the circled-2 to FIG. 5, in a step 50 the controller 14 measures the value of the error signal g. In a step 52 the controller 14 determines whether the absolute value of the error signal g is less than the convergence value calculated in step 36. If the absolute value of the error signal g is found not to be less than the convergence value of the controller 14, the “NO” path is followed from step 52 to a step 54. In step 54, the controller 14 adjusts the control input to the VCO 16 such that the frequency of the signal e₁ is adjusted in proportion to the signal g, but in the opposite polarity to create a negative feedback control system well known to those well-versed in the art.

Steps 50, 52 and 54 are repeated until the absolute value of the error signal g is less than the convergence value. If, during the step 52, the absolute value of the error signal g is determined to be less than the convergence value by the controller 14, the controller 14, in step 56 determines if the closed-loop system has truly converged. For example, if the controller 14 measures the value of the signal h to be less than three times the predetermined value h_(th) then the controller will recalculate the operating parameters and reinitiate the open-loop mode of operation, as depicted by following the circled-3 of FIG. 5 to the step 36 of FIG. 3. If the controller 14 determines that the closed-loop system has properly converged (for example h >3 h_(th)), the controller 14, in a step 58 stops the adjustment of the frequency of signal e₁.

Next, referring to FIG. 6 and FIG. 7 an example of a typical response of the crystal resonant frequency sensor 10 is shown. For this example the crystal arrangement 12 is considered to have a Q of 10,000 and a resonant frequency f_(o) of 10 MHz. The calculated threshold value h_(th) for this example is 200 (the units of FIGS. 6 and 7 are normalized units which represent a digital count). As is shown in the plot of FIG. 6, the crystal arrangement 12 may have several resonant modes. However, once the threshold value h_(th) is achieved the closed-loop mode of operation is initiated and the amplitude of the signal e₂ increases as the output signal e₁ of the voltage controlled oscillator 16 approaches the resonant frequency of crystal arrangement 12.

FIG. 7 shows the frequency error in Hz converging from an initial value of 4 kHz to a final value of approximately zero Hz. The slope discontinuity of the curve in FIG. 7, at about time=1600 microseconds, coincides with the closure of the loop.

Using crystal oscillators to detect the presence of various chemicals is known in the art. For example, in U.S. Pat. No. 5,179,028 (Vali et al.) a high Q crystal oscillator is used to measure concentrations of specific chemicals of interest in the general vicinity of the high Q crystal oscillator. It can be appreciated that the operation of the crystal resonant frequency sensor 10 of this invention can be applied to advantage in the determination of chemical concentrations as described in Vali et al., since the crystal resonant frequency sensor 10 efficiently determines the resonant modes of the high Q crystal oscillator, representing the presence of a chemical of interest. That is, the teaching of this invention can be employed to determine resonator properties which are influenced by an external agency, such as a change in mass due to the presence of a chemical species of interest, a change in temperature, a change in pressure, and other factors that can influence the resonant modes of a resonator.

In FIG. 8, a normalized input signal is plotted with respect to frequency. FIG. 8 illustrates the oscillatory properties of a device under test by a frequency sensor in accordance with the present invention. As shown, there are three points at which the signal goes to zero. The three points represent three, closely spaced resonant frequencies of the device under test. The later two points represent spurious resonance and are generally undesirable. While undesirable, measurements which detail the spurious resonance of a device under test reveal information about the transverse oscillatory properties of the device. Thus, reducing the transverse oscillatory properties improve the design of the device. It is possible, by varying the phase shifter 18, to improve the selectivity of the later two points and thus improve the precision of the estimates of all three frequencies. Such variations might include varying the capacitance of a purely capacitive phase shifter or including other passive elements such as an inductor or a combination of passive elements.

In FIG. 9, a schematic diagram of a two-part mechanical structure whose resonant characteristics can be measured by the present invention is shown. As in the above discussion of the Vali et al. patent, the resonant characteristics of a highly oscillatory system, either mechanical or electrical, can be measured with the crystal resonant frequency sensor 10 and the evaluation method as outlined above. Thus, the teachings of this invention can be employed to determine resonator properties which are influenced by an external agency. That is, the sensor can measure an acceleration a₁ determined by masses m₁ and m₂, springs k₁ and k₂, dampers b₁ and b₂, and an applied force F₁. The applied force F₁, and corresponding acceleration a₁ result in movements x₁ and x₂ of the masses m₁ and m₂, respectively.

Accordingly, an electrical signal from the VCO may be sent to a linear, electromagnetic actuator, such as a BEI model LA-30-27-001 which is capable of producing a desired force for a desired current. The resultant acceleration of the mass could be measured using a variety of accelerometers such as an Entran Model EGCS-A1-2/Z1/L2M/RS, which provides a voltage proportional to acceleration. The output of the phase shifter 12 a may be used to drive the linear actuator. The signal e₂ would be the output of the accelerometer. Otherwise, the operation of the sensor 10 proceeds substantially as described previously.

As an example of the response of this mechanical system, FIG. 10 shows the magnitude of the acceleration of mass m₁ for constant amplitudes of force at various signals. The graph in FIG. 10 showing the response of the mechanical system has more than one resonant peak. Using a purely capacitive phase shifter as before, the expected error signal g may be calculated. This is shown in FIG. 11 for the first peak and FIG. 12 for the second peak. This signal provides a very linear region around the resonant point of the system.

Though the present invention has been described with particular reference to the crystal resonant frequency sensor 10, the teachings of this invention may be used in a precision, linear, crystal-controlled oscillator. In this case, the controller would continue to operate continuously in the closed-loop mode. Here, the convergence value may signal that the circuit was oscillating stably at the desired frequency. The output of the VCO provides the output of the system that could be used for a variety of applications. This approach greatly relaxes the manufacturing tolerances on crystals for use in precision oscillators.

Thus, while the invention has been particularly shown and described with respect to preferred embodiments thereof, it will be understood by those skilled in the art that changes in form and details may be made therein without departing from the scope and spirit of the invention. 

What is claimed is:
 1. A sensor for measuring the frequency of a desired resonant mode of a crystal arrangement, comprising: a voltage-controlled oscillator, having an input and an output, for generating at said oscillator output a signal whose frequency is responsive to a command signal at said input of said oscillator; and a phase shifter having an input and an output, said input of said phase shifter electrically coupled to said output of said voltage-controlled oscillator, said crystal arrangement electrically connected in shunt with said output of said phase shifter; a first multiplier circuit having two inputs, both being coupled with the output of the phase shifter wherein a passive response of said crystal arrangement is multiplied with itself by said first multiplier circuit to generate a first signal having a characteristic amplitude versus frequency curve which peaks substantially at the frequency of the desired resonant mode of the crystal arrangement; a second multiplier circuit having two inputs, one being coupled with the output of the voltage-controlled oscillator and the other being coupled with the output of the phase shifter wherein the two inputs are multiplied by said second multiplier circuit to generate a second signal having a characteristic amplitude versus frequency curve which crosses a predetermined value substantially at the frequency of the desired resonant mode of the crystal arrangement, said crossing occurring within a frequency range in which the first signal is above a predetermined threshold value; and a controller for generating said command signal based on values of said first and second signals; wherein said voltage-controlled oscillator, said phase shifter, said first and second multiplier circuits and said controller form a control loop, wherein said controller adjusts said command signal such that said first signal approaches the predetermined threshold value and said second signal approaches the predetermined value, wherein when said first signal exceeds the predetermined threshold value and said second signal equals substantially the predetermined value the frequency of said output signal of said voltage-controlled oscillator is correlated with the frequency of said desired resonant mode of said crystal arrangement.
 2. A sensor as in claim 1, wherein said predetermined value is zero.
 3. A sensor as in claim 1, wherein said phase shifter is comprised of a capacitor.
 4. A sensor as in claim 1 further comprising: a first filter circuit electrically disposed between said output of said first multiplier and said controller, said first filter circuit being responsive to said first signal for providing a filtered first signal to said controller; and a second filter circuit electrically disposed between said output of said second multiplier and said controller, said second filter circuit being responsive to said second signal for providing a filtered second signal to said controller.
 5. A sensor as in claim 1, wherein said controller is programmable, and further comprising means for programming said programmable controller.
 6. A sensor as in claim 5, wherein said means for programming said programmable controller includes an input device for passing user inputted parameters to said programmable controller.
 7. A sensor as in claim 1 wherein the controller comprises a microprocessor.
 8. A sensor for measuring the resonant frequency of multiple, closely spaced modes of a crystal arrangement, comprising: a voltage-controlled oscillator having an input and an output, said oscillator being adapted to provide an oscillating output signal whose frequency is responsive to a command signal at said input of said oscillator; a programmable phase shifter having an input and an output, wherein an error signal is generated which is proportional to a phase difference between said input of said programmable phase shifter and said output of said programmable phase shifter, said input of said programmable phase shifter electrically coupled to said output of said voltage-controlled oscillator, wherein said crystal arrangement being electrically connected in shunt with said output of said programmable phase shifter; means for programming said programmable phase shifter; a first multiplier circuit having at least two inputs and one output, two of said inputs of said first multiplier being electrically coupled with said output of said programmable phase shifter, wherein a passive response of said crystal arrangement is amplified; a second multiplier circuit having at least two inputs and one output, one of said inputs being electrically coupled with said output of said voltage-controlled oscillator, the second of said inputs being electrically coupled to said output of said programmable phase shifter; a controller having at least two inputs and at least one output, wherein said output is electrically coupled to said input of said voltage-controlled oscillator for providing said command signal to said voltage-controlled oscillator, one of said inputs is electrically coupled to said output of said first multiplier, and a second of said inputs is electrically coupled to said output of said second multiplier; a first filter circuit electrically dispose between said output of said first multiplier and said controller, said first filter circuit having an input and an output and being responsive to said output of said first multiplier applied to said input of said first filter circuit for providing a filtered output signal to said first input of said controller, said output of said first filter circuit defining a signal proportional to an amplitude response of said crystal arrangement to said oscillating output signal of said voltage-controlled oscillator; and a second filter circuit electrically disposed between said output of said second multiplier and said controller, said second filter circuit having an input and an output and being responsive to said output of said second multiplier applied to said input of said second filter circuit for providing a filtered output signal to said second input of said controller, wherein said output of said second multiplier defining an error signal indicative of one of said multiple resonant modes of said crystal arrangement; when in an open-loop mode, said voltage-controlled oscillator, said programmable phase shifter, said first multiplier, said first filter, and said controller form a control loop, wherein said output of said programmable phase shifter is adjusted until said output exceeds a predetermined threshold value; when in a closed-loop mode, said voltage-controlled oscillator, said programmable shifter, said second multiplier, said second filter, and said controller form the control loop, wherein said output of said programmable phase shifter is adjusted until the value of said error signal equals a predetermined value, whereby the frequency of said oscillating output of said voltage-controlled oscillator is correlated with the frequency of one of said multiple modes of said crystal arrangement.
 9. A method for determining the frequency of a desired resonant mode of a crystal arrangement, comprising the steps of: calculating operating parameters related to the desired resonant mode of the crystal arrangement, the operating parameters comprising a frequency range of operation, a final error value, and a threshold value; providing a control signal to a voltage-controlled oscillator, wherein the frequency of a sinusoidal output signal of the voltage-controlled oscillator is responsive to the control signal; and providing the sinusoidal output signal of the voltage-controlled oscillator to the crystal arrangement wherein the frequency of the sinusoidal signal being within the calculated frequency range, the sinusoidal output signal first passing through a phase shifter having an input and output, wherein an error signal is generated which is proportional to a phase difference between the input of the phase shifter and the output of the phase shifter, when in an open-loop mode, measuring a response of the crystal arrangement to the sinusoidal output signal; in the open-loop mode, adjusting the frequency of the sinusoidal output signal until the value of the measured response of the crystal arrangement is greater than the calculated threshold value; when in a closed-loop mode, measuring the error signal; and in the closed-loop mode, adjusting the frequency of the sinusoidal output signal until the value of the measured error signal is equal to the calculated final error signal value, whereby the frequency of the sinusoidal output signal is the frequency of the desired resonant mode of the crystal arrangement.
 10. A method as in claim 9, wherein the step of calculating operating parameters further comprises the step of calculating a system gain, wherein the rate at which the frequency of the sinusoidal signal can be adjusted is increased such that less time is required to determine the frequency of the desired mode of the crystal arrangement.
 11. A method as in claim 9, wherein the step of adjusting the frequency of the sinusoidal signal until the value of the measured response of the crystal arrangement is greater than the calculated threshold value further comprises the step of adjusting the control signal to the voltage-controlled oscillator, wherein the frequency of the sinusoidal output of the voltage-controlled oscillator sweeps within the calculated frequency range, and wherein the value of the response of the crystal arrangement approaches the calculated threshold value.
 12. A method as in claim 11, wherein the step of adjusting the frequency of the sinusoidal signal until the value of the measured error signal is equal to the calculated final error signal value further comprises the step of adjusting the control signal to the voltage-controlled oscillator wherein the frequency of the sinusoidal output of the voltage-controlled oscillator is finely adjusted within the calculated frequency range, and wherein the value of the response of the crystal arrangement approaches the calculated final error value.
 13. A method for determining the resonant frequency of multiple, closely spaced modes of a crystal arrangement, comprising the steps of: calculating operating parameters related to one of the multiple, closely spaced resonant modes of the crystal arrangement, the operating parameters comprising a frequency range of operation, a final error value, and a threshold value; providing a control signal to a voltage-controlled oscillator, wherein the frequency of a sinusoidal output of the voltage-controlled oscillator is responsive to the control signal, and the frequency of the sinusoidal output of the voltage-controlled oscillator is within the calculated frequency range; providing the sinusoidal output of the voltage-controlled oscillator to the crystal arrangement, the sinusoidal output first passing through a phase shifter having an input and an output, wherein an error signal is generated which is proportional to the phase difference of said input of said phase shifter and said output of said phase shifter; when in an open-loop mode, measuring a response of the crystal arrangement to the sinusoidal output; in the open-loop mode, adjusting the control signal to the voltage-controlled oscillator, wherein the frequency of the sinusoidal output of the voltage-controlled oscillator sweeps within the calculated frequency range, and wherein the value of the response of the crystal arrangement approaches and exceeds the calculated threshold value; when in a closed-loop mode, measuring the error signal; and in the closed-loop mode, adjusting the control signal to the voltage-controlled oscillator wherein the frequency of the sinusoidal output of the voltage-controlled oscillator is finely adjusted within the calculated frequency range until the value of the measured error signal is equal to the calculated final error signal value, whereby the frequency of the sinusoidal output is the frequency of one of the multiple, closely spaced resonant modes of the crystal arrangement. 